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Any nonlinear equation of the form y" = 'Ylm=o ^n{z)y"' has a (generally 
branched) solution with leading order behaviour proportional to {z — zq)~'^/^^~^'^ 
about a point zq; where the coefficients a„ are analytic at zq and ajv(zo) 7^ 0. We 
consider the subclass of equations for which each possible leading order term of this 
form corresponds to a one-parameter family of solutions represented near zq by a 
Laurent series in fractional powers oi z — zq. For this class of equations we show that 
the only movable singularities that can be reached by analytic continuation along 
finite-length curves are of the algebraic type just described. This work generalizes 
previous results of S. Shimomura. The only other possible kind of movable singular- 
ity that might occur is an accumulation point of algebraic singularities that can be 
reached by analytic continuation along infinitely long paths ending at a finite point 
in the complex plane. This behaviour cannot occur for constant coefficient equa- 
tions in the class considered. However, an example of R. A. Smith shows that such 
singularities do occur in solutions of a simple autonomous second-order differential 
equation outside the class we consider here. 



1 Introduction 

The most fundamental result in the theory of differential equations in the complex domain 
is Cauchy's existence and uniqueness theorem. 

Theorem 1 ( Cauchy) Let /i, . . . , /„ he analytic functions in a neighbourhood of the point 
(2:0, r^i, . . . , r/n) in C"^^. There is a unique n-tuple of functions yi,...,yn analytic in a 
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neighbourhood of zq ^ C such that for all z G ^l, 

y[ = fiiz; yi,...,yn) and y,{zo) = rji, for all i = 1, . . . ,n. 

Cauchy's theorem is a local result. It guarantees the existence of a local solution 
of a differential equation provided that the equation is well behaved (analytic) at the 
initial condition. However, it says nothing about the nature of this solution after analytic 
continuation. For example, Cauchy's theorem guarantees that the initial value problem 

has a unique solution in a neighbourhood of 2; = 0. This solution is 

y{z)=c[{l + z)/il + c'z)Y/'. 

For c 7^ this solution is analytic on {z : \z\ < min{l, |c|~^} but has branch points 
at z = — and z = —1. The location of the singularity aX z = — varies with the 
initial condition. Such singularities are called movable. The location of the singularity 
at 2; = — 1 does not depend on the initial condition. Such singularities are called fixed. 
Heuristically speaking, fixed singularities occur at points where the equation itself is sin- 
gular in some sense, whereas movable singulaities develop spontaneously. All singularities 
of autonomous equations are movable. 

A singularity of a solution of the linear differential equation 

+ a„„i(^)i/("-i) + ■ ■ ■ + a^{z)y' + a,{z)y = f{z) 

can only occur at a singularity of one of the coefficient functions oq, . . . , a„_i or /. It 
follows that linear equations do not have movable singularities. In [22], Painleve showed 
that the only movable singularities of solutions of first-order ODEs of the form 

F{z;y,y')=0, (1) 

where F is a polynomial in y and y' with coefficients that are analytic in some common 
domain, are poles or algebraic branch points. The situation for higher-order equations is 
much more complicated and very poorly understood. For example, the general solution 
of the equation 

iyy" - y"? + ^yy" = o 

is y{z) = cexp{{z — 2:0)^^}, which has a movable essential singularity at zq. Painleve 
considered the equation 

y" = %i^y", (2) 

which has the general solution 

y{z) = tan(log(ci2; - C2)), (3) 
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where ci and C2 are constants. If ci 7^ 0, this solution has movable poles accumulating 
at the movable branched singularity at z = c^jcx. Painleve's example is often used 
as a warning about the overzealous use of so-called Painleve tests. Using formal series 
methods, it is very easy to show that equation ([2]) has a two-parameter family of Laurent 
series solutions about poles. Although every nontrivial solution contains infinitely many 
movable poles, it also contains a branched accumulation point of poles, which is not 
detected by the standard analysis. 

Of particular relevance to our investigation is the example of Smith [29] , 

l/" + %V + Z/ = 0. (4) 

He showed that if a solution of (jlj) can be continued analytically along a curve 7 of finite 
length up to but not including the point zq? then y has an algebraic branch point at 
Zq. Moreover, he showed that there are solutions such that analytic continuation along 
a curve of infinite length ends at a point in the finite 2;-plane that is an accumulation 
point of these algebraic branch points. This kind of accumulation of singularities is quite 
different from that in Painleve's example. We will discuss this example in more detail in 
section HI 

For third order equations other phenomena arise. For example, the general solution 
of the Chazy equation 

III n II o '2 

V =^yy - 3y 

has a movable natural barrier [H |5| . 

In this paper we will primarily be concerned with movable singularities of solutions of 
the equation 

N 

/(^) = ^a„(%(^r, nGN, n>2. (5) 

n=0 

A singularity at zq of a solution of equation ([5]) is called movable if ao, . . . , oat are analytic 
at zq and aN^zo) 7^ 0. Note that the case < < 1 is linear, in which case y will have 
no movable singularities. For the rest of this paper we assume that N > 2. The main 
result of this paper is the following theorem, which is a characterization of the movable 
singularities of equation It is a generalization of results of Shimomura described 
below. 



Theorem 2 For N > 2, suppose that there is a domain Q C C such that oq, . . . , ctAr are 
analytic and that aN{zo) ^ on Q. Suppose further that for each zq ^ Q and for each cq 
such that 

9 ^r-L^ 

(6) 



_ 2 N + 1 
'° aNizo)iN-ir 



equation ([5]) admits a formal series solution of the form 



j=0 



Then 
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1. 



For each cq satisfying ([6]) and for each f3 E C, there is a unique formal series 
solution of the form ((Tj) such that C2(Ar+i) = (3. 



2. Given cq and C2(Ar+i) as above, the series (JTj) converges in a neighbourhood of zq. 

3. Now let y be a solution of equation ([5]) that can be continued analytically along a 
curve 7 up to but not including the endpoint zq, where the coefficients aj are analytic 
on 7 U {2:0} and is nowhere zero on j U {zq}. If •y is of finite length, then y has 
a convergent series expansion about Zq of the form ([7]). 

4- If y cannot be represented by a series expansion about zq of the form (I7l) then 7 is 
of infinite length and zq is an accumulation point of such algebraic singularities. 

By "accumulation point" in part 4 of the theorem we mean that given any e > there 
exists a straight hne segment I in the disk of radius e centred at zq with endpoints -Zi G 7 
and Z2 such that analytic continuation of y along 7 up to zi and then along I ends in an 
algebraic singularity at Z2. 

We will see in section 2.3 that the assumptions of theorem [2] correspond to one or two 
differential relations relating the coefficients of equation ([5]) depending on whether N is 
even or odd respectively. These relations are resonance conditions. That is, they represent 
the vanishing of an obstruction that occurs when the recurrence relation determining the 
coefficients in the series ([7]) breaks down. We will see in section 2.1 that equation (jS]) can 
be normalized so that = 2{N+1)/{N — 1)'^ and a^-i = 0. With this normalization the 
resonance condition in the even N case is a%_2 = 0. In the odd N case the two resonance 
conditions are equivalent to a^_2 = and another condition, which is listed below for the 
first few odd values of N: 



The fact that the resonance conditions are satisfied if each a„ is a constant is apparent 
from the following "explicit" integration. If y is non-constant then multiplying equation 
(IS]) by 2y' and integrating gives 



for some constant k. Equation ([8]), which can be solved in terms of hyperelliptic functions 
or their degenerations, is of the form ([T]). Therefore, we know that the only singularities of 
its solutions are algebraic, so the accumulation of singularities described by part 4 of the 
theorem does not occur here. It remains an open question whether such an accumulation 
of singularities can occur for equations of the form iQ. This phenomenon is known to 
occur in the example of Smith (jl]), which lies outside of the class considered here. 

The role of movable singularities is particularly important in the theory of integrable 
systems in the context of the Painleve property. An ODE is said to possess the Painleve 



N = 3: a'oiz) = 0, 
A^ = 5: [4ai{z) - al{z)]' = 0, 
N = 7: [10a2(z) - 904(2)05(2)]' = 0. 




n=0 
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property if all movable singularities of all solutions are pole^. If solutions of an equation 
possess no fixed singularities in the finite plane, e.g., if the equation is autonomous, then 
the Painleve property is equivalent to the property that all solutions are meromorphic. 
The connection between the integrability of an ODE and the singularity structure of its 
solutions appears to have been first exploited by Kowalevskaya, who used it to find a new 
integrable case of the equations of motion of a spinning top [HI [17] . 
In [8], L. Fuchs studied equations of the form 

y' = Fiz;y), 

where F is rational in y and with coefficients that are analytic in some domain. He showed 
that the only nonlinear equation in this class with the Painleve property is the Riccati 
equation, which corresponds to the case in which F is quadratic in y. The general solution 
of the Riccati equation 

y' = a{z)y^ + b{z)y + c{z), (9) 



where a ^ 0, is given by 



y{z) = -^^, 10 



where w is the general solution of the linear equation 

a(z)—— — W(z) + a(z)b(z)] — h c(z)a'^(z)w = 0. 

dz^ dz 

Since the singularities of w are fixed, we see from the transformation fllOl) that the only 
movable singularities of y are poles, corresponding to zeros of w. Here we see our first 
example of simple movable singularities (poles) detecting integrable (in this case, lineariz- 
able) equations. 

Painleve, Gambler and Fuchs classified all equations with the Painleve property of the 
form 

y" = F{z;y,y'), 

where F is rational in y and y' with coefficients that are analytic in some domain. They 
found that each such equation could be transformed to one of fifty canonical forms. Six 
of these canonical equations are now called the Painleve equations Pj-Pvi- All of the 
other canonical equations could either be solved in terms of classically known functions, 
by quadrature, by solving linear ODEs or in terms of solutions of one of Pi-Pyi- The 
Painleve equations were subsequently shown to be compatibility conditions for certain 
isomonodromy (linear) problems. They are extremely important equations in the theory 
of integrable systems. The solutions of the Painleve equations are often referred to as 
nonlinear special functions. For each Painleve equation, all solutions are meromorphic on 
the universal cover of C \ 5*, where 5 is a finite (possibly empty) subset of C. 

Interest in the Painleve property and movable singularities was rekindled in the late 
twentieth century through connections with the spin-spin correlation function of the two- 
dimensional Ising model (Wu, McCoy, Tracy and Barouch ^ST[ fTU]). holonomic quantum 



^Some authors say that an ODE is said to possess the Painleve property if all solutions are single- valued 
around all movable singularities. 
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fields (Sato, Miwa and Jimbo [25]) and the theory of sohton equations. In their studies 
of the asymptotic behaviour of solutions of soliton equations Ablowitz and Segur [3] , and 
Ablowitz, Ramani, and Segur [1], '2] discovered that ODE reductions of soliton equations 
often led to equations of Painleve type, possibly after a transformation of variables. 

Relatively simple tests of strong necessary conditions for the Painleve property can be 
used to detect integrability [HIS]- These tests involve showing that formal Laurent series 
solutions exist and that certain resonance conditions are satisfied. However, it is a much 
more difficult task to show that a given equation actually possesses the Painleve property, 
without having a more-or-less explicit representation of the solution. 

The main idea of many of the standard proofs of the Painleve property is to character- 
ize the singular Laurent series expansions in terms of locally analytic variables satisfying 
a regular initial value problem. The main tool in this respect is the following result, which 
can be found in a number of sources [lOl [11] and follows from the fact that Cauchy's ex- 
istence and uniqueness theorem can be strengthened to give a lower bound on the radius 
of convergence of the solution. 

Lemma 3 (Painleve) Let /i, . . . , /m be analytic functions in a neighbourhood of the point 
{a, rji, . . . , rjm) in C™^"^. Let 'j be a curve with end point a and suppose that yi is analytic 
on j\ {a} for i = 1, . . . ,m and satisfies 

y'i = fi{z;yu.. .,ym)- 

Let (zn) be a sequence of points such that Zn & 'j, Zn ^ a and yi{zn) rji as n ^ oo , 
for all i = 1, . . . ,n. Then each yi is analytic at a. 

Painleve himself provided a proof that the first Painleve equation 

y" = Qy^ + z (11) 

possesses the Painleve property. This proof, which appears in a number of forms in the 
literature (e.g., Ince [H], Golubev [9]), had a number of gaps in it that have been filled 
by several authors (Hukuhara [131 [21], Hinkkanen and Laine [12], Shimomura [26] )• AH 
of these proofs proceed in the following manner. 

i. Show that if a solution y of equation (fTTl) has a pole aX z = Zq then it is a double 
pole and the coefficient j3 oi {z — zqY is arbitrary. Fixing a value for j3 uniquely 
determines all the other coefficients. 

ii. Consider a finite length curve on which y is analytic except at the endpoint Zq where 
it is singular. Show that y is unbounded on 7. Show that if y is not bounded away 
from zero on 7 then 7 can be deformed to a new finite length curve 7 ending at zq 
so that y is bounded away from zero on 7. 

Hi. Using the form of the series expansion in introduce a new variable u := y"^''^, 
for some choice of branch and another variable v such that v{zq) encodes the value 
of the resonance parameter (3. The functions u and v satisfy a regular initial value 
problem with u{zq) = and v{zq) = k for k G C. 
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iv. A function W of z, y and y' is shown to be bounded on 7. This function is in some 
sense an approximate first integral of the equation. 



V. If A := hminf^9^_^oo \y{z)\ is finite and positive then the boundedness of W leads 
to Z/i = 2/ and 7/2 = u' satisfying the conditions of lemma [3l so y is analytic at Zq. 

vi. li A = 00 (i.e. if lim^g^^^oo = 00), then the boundedness of W (and the 

correct choice of branch in the definition of u) shows that v must be bounded on 7, 
so applying lemma [3] to the initial value problem for u and v mentioned in in shows 
that y must have a double pole at 2:0- 

Many of the other methods for proving that the Painleve equations possess the Painleve 
property rely on the underlying isomonodromy problems and Riemann-Hilbert techniques 
(Miwa [20], Malgrange [18], see also Fokas, Its, Kapaev and Novokshenov [7]). These 
methods explicitly exploit structure that is closely associated with the integrability of the 
Painleve equations and they are therefore unlikely to generalize directly to the noninte- 
grable equations that we consider. Other proofs of the Painleve property for the Painleve 
equations are Steinmetz [30], who uses differential inequalities, Erugin [6], and Joshi and 
Kruskal [15]. 

Shimomura [271 [28] studied the quasi- Painleve property, which is also known as the 
weak Painleve property (see [231 121] )• An ODE is said to possess the quasi-Painleve 
property if all movable singularities of all solutions are at most algebraic branch points. 
He proved that for any G N, analytic continuation of any solution y of 



along a finite length curve ends at a point where either y is analytic, has a pole or has an 
algebraic branch point. The method of proof is essentially the same as outlined above for 
equation (II ip except that the initial value problem for the new variables u and v described 
in in is now for z{u) and v{u), rather than u{z) and v{z). In this way it is shown that z 
is an analytic function of u but its inverse is in general analytic in a fractional power of 

Z - Zq. 

The main result of the present paper is a generalization of Shimomura's result. Due to 
the complexity and generality of the series expansions we would otherwise have to consider, 
we show directly that the existence of algebraic formal series solutions is equivalent to 
the existence of a bounded function W (modulo the same curve modification arguments 
outlined in the description of the proof of the Painleve property in ii). The choice of 
suitable functions u and v comes directly from the expression for W . The equations 
for W , u and v are implicit in the sense that an algorithm is presented for calculating 
their coefficients in terms of the functions Oj and that these algorithms are well defined 
provided that the resonance conditions are satisfied. This is especially useful in the odd 
N case in which one of the resonance conditions is determined recursively. It is this use 
of representing certain associated equations using recurrence relations that allows us to 
get a results for such a large class of equations. The rest of the proof follows the general 
pattern above. 



y" = 




(12) 
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The possibility that algebraic branch points could accumulate along infinite length 
curves in a bounded region of the complex plane appeared in Smith's study of equations 
of the form 

y" + fiy)y' + 9iy) = P{z), 

where / and g are polynomials and the degree of / is greater than the degree of g and P 
is analytic in some domain [29]. 

In section 2 we will transform equation ([5]) to a canonical form that will simplify future 
calculations. We will then show how to construct singular series expansions of solutions. 
The different nature of the leading orders and the resonance conditions in the even and 
odd cases will be described. In section 3 we show that the existence of these singular 
series expansions implies the boundedness of a certain function W. This section then 
follows closely the outline above. In section 4 we consider analytic continuation along 
infinite length curves in a bounded region. We end with a discussion in section 5. 



2 Formal series expansions 

In this section we use standard methods to determine when formal Laurent series ex- 
pansions of solutions of equation ([5]) in fractional powers of z — zo exist. The method 
described in this section is the most direct way of verifying whether a given equation of 
the form ([5]) satisfies the assumptions of theorem [21 The fundamental difference between 
the even N and odd cases will be immediately apparent. Although this section only 
deals with formal series expansions, it will be shown at the end of section 3 that they 
converge. 



2.1 Canonical form of equation ([5]) 

Suppose that qq, ■ ■ ■ ,aN are analytic in a neighbourhood of zq and that aN{zo) ^ 0. We 
begin by transforming equation ([5]) to a canonical form in order to simplify our analysis. 
To this end, let 



{N-iyaN{z) J ' / Najviz) 

yiz) = f{z)y{~z) + g{z), ~z = f f-^r) dr. 



This transformation is analytic and invertible in a neighbourhood oi z = zq. Under this 
transformation, z = zq is mapped to z = and equation ([5|) is mapped to an equation 
of the same form with z and y{z) replaced by z and y{z) and with aj{z) replaced by 
aj{z), j = 0, . . . , N, where all the coefficients dj are analytic in a neighbourhood Cl of 0, 
ajv(5) = 2(A^ + 1)/(A^ — 1)^ and aiv_i(5) = 0. Throughout sections 2 and 3 we will replace 
equation by this canonical form. By a slight abuse of notation, we will not use tildes 
on the new variables but replace them by the orginals. Hence, without loss of generality. 
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we restrict our attention to the equation 

y"{z) = ^n{z)y{z) + -A^^,y\^). (13) 

n=0 ^ ' 

in an open set VL in which the coefficients a„ are analytic. 
2.2 Leading order analysis 

Choose a G fi. We will look for formal solutions of equation f|T3l) of the form 

oo 

y{z) = Y,Cr.r-', (14) 
n=0 

where p and q are positive, Cq 7^ and C, = z — a. Substituting the expansion (IT^ into 
equation f|T3|) and keeping only the leading-order terms (i.e., the terms that are largest 
for small Q on both sides, we have 

g(g + i)coC--^ + ■ ■ ■ = ^T^^^C-^ + ■ ■ ■ , C - 0. 

Equating powers of Q in the leading order terms gives q = 2/{N — 1). Equating the 
coefficients of (^-2A^/(a^-i) then gives Cq~^ = 1. So the leading order behaviour of y is 

y^CoC'^(^-'\ cr' = l. 

In terms of the original variables in equation Q, this condition on cq corresponds to 
equation ([6]). 

It is at this point that we ffist see that the parity of plays an important role. If 
is even then the mapping ( 1—^ ^-2/{Ar-i) j-^^^g _ ^ branches about C = 0- The different 
choices of Cq correspond to different choices of branch. The choice of Cq can be effectively 
absorbed into the choice of branch, so we can take Co = 1, i.e., y{z) ~ (-2/(a^-i). if 

= + 1 is odd then the mapping C ^ (-^/(^-i) = ('^^ has K = {N - l)/2 
branches, whereas there are — 1 = 2K choices of cq. For some e > 0, let xq be one 
of the possible values of the function f{z) = z^^l^ at z = e. Analytic continuation of 
/ in a clockwise direction n times around the circle |z| = e gives the value uj'^Xq, where 
uo = exp(27ri/i^) is the primitive i^th root of unity. The choices for cq naturally fall 
into two classes: those for which Cq = 1 and those for which Cq = —1. For a fixed 
determination of , it can be seen that analytic continuation around C = cannot 
change the effective value of Cq from one class to the other because analytic continuation 
effectively only multiplies Cq by a power of u. Hence, absorbing the choice of Cq as much as 
possible into the choice of branch, we see that the odd A^ case has two essentially different 
possible leading order behaviours. 
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2.3 Resonance conditions 



Having computed the possible leading order behaviours, we now consider how they can be 
extended to series solutions. A natural guess would be to extend them as Laurent series 
in (^2/(Af~i)^ jj^ j^Yiis case the left side of equation (|T3l) would also be a Laurent series in 
^2/(7v-i)_ However, ( can be expanded as a Laurent series in (^/(^-^) about if and only 
if is odd (in which case ( = ((^2/(Af-i)'j(Af-i)/2'j_ Therefore, the right side of equation 
f fT3|) would also be a Laurent series in is odd, but not if is even and some 

of the coefficient functions oq, . . . , aN-2 are not constant. However, if A^ is even and we 
extend the leading order term to a Laurent series in then both sides of equation 

( IT3|) will also be Laurent series in this variable. 

Proof of theorem [2] part 1 

We look for formal series solutions of equation f|T3|) of the form 



viz) =J2^nC''-', (15) 
n=0 

which corresponds to (fT4l) with p = 1/{N — 1) and q = 2/{N — 1). For even A^ we have 
Co = 1. For odd A^ we have Cq^"^'*'^^ = ±1 and C2k+i = for all = 0, 1, 2, . . . Substituting 
( IT5|) into equation ( |T3l) and equating coefficients of like powers of ( gives a recurrence 
relation for the coefficients c„. For any positive integer r, the lowest power of ( such that 
Cr appears as a coefficient on either side of equation (fT3|l is := — 2N)/{N — 1). The 
coefficient of on the left side of f|T3|) is 

{r -2){r - N - 1) 
(AT -1)2 

while the corresponding coefficient on the right side has the form 

2N{N + 1) 



'Cr ~\~ -Pj.(Co, . . . , Cf— l). 



where Pr is a polynomial in its arguments. Equating these coefficients gives a recurrence 
relation of the form 

(r + A^ - l)(r - 2Ar - 2)c, = (A^ - l)2p,(co, . . . , c,_i), (16) 

for all r = 1, 2, . . . 

Equation ( |T6l) is a recurrence relation for the coefficients Cr with a resonance at 
r = 2(A^ +1)- In other words, we can use ( |T6i) to determine inductively in terms 
of Co, . . . ,Cr-i except when r = 2{N + 1). A necessary and sufficient condition for the 
existence of a formal series expansion of the form (ITS!) is that P2(n+i){co, ■ ■ ■ , Cr-i) = 0. 
If this condition is satisfied then C2(Ar+i) is arbitrary. □ 

In the even A^ case there is effectively just one leading order behaviour, so we get just 
one resonance condition. We will see later that this condition is cl'^^2{'^) = 0- Since a E Q 
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is arbitrary, this becomes the differential condition ci^_2 = 0. That is aN-2{z) = Az + B 
for some constants A and B. In the odd N case, there are two leading order behaviours 
leading to two resonance conditions. These two conditions taken together lead to the 
condition a'l^_2 = and another condition, which we express in a more convenient way 
(in equations and which are still not in closed form) in the next section. 



3 Movable algebraic branch points 

We begin this section by establishing the properties of certain functions that will be useful 
in the proof of part 3 of theorem [2l Part 2 of theorem [2] will be a consequence of the 
existence of a particular regular initial value problem that arises in the proof of part 3 of 
theorem [21 



3.1 A bounded function 

Let y be a solution of equation f|T3l) that can be analytically continued along a finite 
length curve 7 up to but not including the endpoint zq. The purpose of this section is to 
identify a function W{z; y, y') that is bounded along 7 \ {2:0} whenever y is bounded away 
from on 7 \ {zq}. If the coefficients a„ are all constants then a suitable function would 
be the first integral 

N+l 



k 

k=l 



More generally, we consider the ansatz 



, / N \ N+l 



Wiz) := y'(zf + ( E ^ j y'(z) - 2 E ^/(^). (17) 



k=l 



where the functions &i, . . . , &Ar-i are analytic at zq but yet to be determined. 

We will now obtain an integral representation for W . Differentiating equation ( JT71) 
with respect to z and substituting the second derivative of y{z) from equation f|T3|) and 
y'"^ from f|T7j) . we get a linear equation of the form 



W + P{z, l/y)W = Q{z, l/y)y' + R{z, l/y) + S{z, y), 

where 

7V-1 

P(z, l/y) = J2 kbk{z)/y'^\ 

k=l 

N-l /N-l \ /N-l 

Q{z, l/y) = b',{z)/y' + kbk{z)/y'^' Y ^^^^1^" 

k=l \k=l / \k=l 

and R{z,l/y) + S{z,y)= (19) 

N-l N N+l /N-l \ /N+l \ 

Y bk{z)/y' Y «^ (^)^' - 2 E <-ii^)y'lk - 2 E khu{z)/y'^^' ^ «^-i^V^ • 

k=l j=0 k=l \k=l / \k=l J 
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The functions P, Q, R, S are polynomials in their second arguments. We separate non- 
positive powers of y{z) into R{z, 1/y) and positive into S{z,y). Equating the coefficients 
of the positive powers of y{z) in (fTOll gives 



N-l , ,N-n _ I . \ 2 \ 

s{^^y) = E ( ( E ^ + ^ + (20) 

n=l ^ ^ 171=1 ^ ' 

Now consider the case in which N is even and use the recurrence relation 

m=l 

to define uniquely &i, . . . , &Ar_i, so that S{z, y) = 0. Note that the coefficient of hn{z) on 
the left side of equation (I2T1) does not vanish for n = 1, . . . , — 1 since is even. 

For any point z on 7, we denote by 7^ the part of 7 up to 2; (in particular, 7 = 720). 
Using the integrating factor 



E{z) := exp (^J P(C,lMC))dc), (22) 



"Yz 

the solution of equation (|T8l) can be written as 

W^(^) = ^['' + 1 {^(^' Vy(r))y' + ^(r, l/t/(r))} i5;(r) dr^ , (23) 
for some constant k. 

We now impose the main assumption of theorem [2], namely that there is a formal 
series solution y of equation f|T3|) that is a Laurent series in (z — zoY^^'^~^^ with leading 
order behaviour given by ?/ ~ (z — zq)"^''^^""^-'- For such a solution, -^(-z), (5(-2, l/y{z)) and 
i?(2;, l/?/(z)) are power series in {z — zoY^^^~^\ Furthermore, E{z) = 1 + (^{z — zq) n-i^ 

and Q{z, l/y{z))y' = — ^^^"^ {z — zq)^^ + • ■ ■ . Hence the integral in equation fl23l) shows 
that the series expansion for W about zq contains a logarithm if b'i{zo) 7^ 0. However, from 
the definition f|T7j) . W has a Laurent series expansion in {z — zqY^^^'^K So b[{zo) = 0. 
Furthermore, this analysis can be repeated for any point zq G Q, therefore b'l = 0. On 
substituting n = 1 in the recurrence relation ( 12T1) . we see that bi = a'j^_2- So we have 
proved the necessity of the condition a'^_2 = 0. 

Next we consider the case in which N = 2K + 1 is odd. We define the first 
K = {N — l)/2 functions bi,...,bK using the recurrence relation (I2T]) . From equation 
( I20I) . we see that S{z,y) is a polynomial in y of degree at most K = {N — l)/2. The 
coefficient of y^ in S{z, y) is 



^ m=l ' 



By assumption there are two families of solutions at zq that are Laurent series in [z—z^yf^ 
with asymptotic behaviour y ~ cq{z — zq)"^^'^- The two families correspond to the 
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cases Cq = 1 and Cq = —1. For such solutions, P, R and E are all power series in 
{z—zqY^^, where E is again an integrating factor given by (122|) . We also have the following 
leading-order behaviours: S{z,y{z)) ~ pcQ {z — Zq)~^ , Q{z, l/y{z))y' ~ —^^^^^{z — Zq)~^ , 
E{z) ~ 1. Integrating equation (ITSl) we have 

W{z) = ^(k+ [ {Q{t, l/yirW + R{r, l/vir)) + S{t, y{r))} E{r) dr) . (24) 



7z 



E{z 

The coefficient of (2; — 2^0)"^ in the integrand is now 

1{zo]Cq) := p{zo)cq 



iV- 1 ' 



Hence to avoid a logarithm in the series expansion of W in both the 1 and 

Cq = —1, we must impose the conditions l{zo] 1) = l{zo; —1) = 0. These conditions are 
equivalent to b[{zo) = and p{zq) = and must hold for all zq in a nonempty open set, 
leading to a%_2 = and 

^a'j^-i{z) -IY] ! I ^ bm{z)aK+Uz) = 0. (25) 
K 2 ^-^ A + 1 + m 

m=l 

Furthermore, knowing that equation ([25D is satisfied, we can now choose hx+x arbi- 
trarily (e.g., hK^i = 0) and use the recurrence relation to define the functions bk for 
k = K + 2, . . . , N — 1, which means that S{z, y) = 0. 

We have shown in both the even and odd cases that W defined by ffTTj) . where 
the functions bk are defined by the recurrence relation (12T!) and b[ = a'^_2 = 0, is given 
by equation (!23|) for some constant k. Now let 7 be any curve in Q of finite length L 
ending at 2;o ^ ^ and let y be any function analytic on 7 \ {^o} such that there exists a 
constant B > 1 and |y(2;)|^^ < B for all z G 'y\{zQ}. Note that we are no longer assuming 
anything about any possible singularity of y ai zq. 

Since 7 and its endpoints are in Q, there is a constant A and functions p2, ■ ■ ■ ,Pn] 
12, ■ ■ ■ , (l2N-i]ro, . . . , riy_i that are all analytic and bounded by A on 7 such that 

N 2N-1 N-l 

P{z, 1/y) = Y.pu/y\ Q{z, i/y) = ^kly' R{z, l/y) = n/y". 

k=2 k=2 k=Q 

So on 7, 

N N 

\p{z, i/y{z))\ < J2 \Pk\/\y\' <J2ab'<{n- 1)AS^ 

k=2 k=2 

and similarly 

\R[z,l/y{z))\ < NAB^~\ 
So from equation (1221) and the inequalities exp(— |x|) < | expx| < exp(|x|), we have 

exp(-(A^ - 1)AB^L) < \E{z)\ < exp((iV - 1)AB^L). (26) 
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Therefore it follows from equation (l23l) that in order to prove that W is bounded on 7, it 
is sufficient to show that 



I := [ Q{T,l/y{r))y'E{T)dT 



is bounded. Now 



= -E_/ 



fc=2 
2N-1 



A:=2 
2N-1 



, V- / gU^) + gfc(^)nr,l/l/(r)) E{t) 

for some constant C. The integrand is bounded on 7, therefore I is bounded and hence 
so is W. 

In summary, we have proved the following. 

Lemma 4 For N > 2, suppose that there is a domain Q G C such that oq, . . . , aN-2 o'^e 
analytic on il. Suppose further that for each zq E and for each cq such that Cq~^ = 1 
equation f fT3|) admits a formal series solution of the form 

00 

j=0 

Then 

a"N-2 = (27) 

and there are functions (b^)^^^ that are analytic on Vl and satisfy the recurrence relation 
( !2T|) for n = 1, . . . , N — 1. Furthermore, let 'f be a finite-length curve in Q and let y be a 
solution of equation ( IT3l) that is analytic and bounded away from on 7, then the function 
W defined by equation ( IT71) is bounded on 7. 



If is even then the recurrence relation (1211) defines the functions 6^ uniquely with- 
out further constraints. However, if N is odd then the left side of (pTj) vanishes when 
n = {N + l)/2. This gives the second constraint (!25|) . If this constraint is satisfied, then 
the function 6(7v+i)/2 is not determined by the recurrence relation (pT!) and can be chosen 
to have any value (e.g., zero). In section [231 we saw that the existence of formal Laurent 
series expansions in fractional powers oi z — zq was equivalent to one or two resonance 
conditions being satisfied in the even and odd cases respectively. Lemma H] shows 
that equation (127]) (in the even and odd cases) and equation fl25|) (in the odd A^ case) 
are necessary consequences of these resonance conditions. We will see subsequently that 
these conditions are also sufficient. 
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3.2 Proof of theorem [2] part 3 

Let yhe a. solution of equation (fT3l) that can be continued analytically along a finite length 
curve 7 up to but not including the endpoint zq, where the coefficients aj are analytic on 
7 U {2:0}. If y is bounded on 7 then so is 

N 



y' = c+ [ X^«n(cMcrdc. 

•Jy n 



"y n=0 

Hence, from lemma [HI y is analytic at zo, which is a contradiction. So we have shown that 
y is unbounded on 7. In other words, limsup^g^^^o \y{^)\ = c>o. Following the main ideas 
behind Painleve's work on the Painleve property, we divide the proof into three parts 
depending on the value of 

A := liminf 

Case 1: < A < 00. 

Since y is bounded away from on 7 near zq, lemma H] shows that W is bounded. By 
assumption, there exists a sequence of points z„ G 7 such that Zn — ^ Zq and y{zn) — > r] 
for some i] E C. Viewing equation f|T71) as a quadratic equation in y', we see that there is 
a subsequence of (zn) on which y' also approaches a finite value. Hence, applying lemma 
|3] to equation f|T3|) with yi = y and 2/2 = y', we see that y is regular at zq, which is a 
contradiction. 



Case 2: A = 00. 

In this case ?/ — > cxd as 2; 

for y' gives 



1 



2/ = -1 




on 7 and W is bounded on 7 near zq. Solving equation ([T 

1/2 



E 



A:=l 



A; 



, (28) 



for some choice of branch of the square root. 

Assume for now that is even. We introduce a new variable u such that y = l/u"^. 
There is a choice of branch if we take this as a "definition" of u, which we will fix shortly. 
Define the function F{z; u), which is analytic for u near zero and z near zo, by 

2 



F(z-u) 



u 



2(N+l) 



N+l 



2k 



+«E 



. fc=i 



k=l 
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N 



{N 



Qfc-l(^) 2(N-k+l) 

Uj 



Qfc-l(^) 

N-l 



k=l 



k 



+ bk{z)bi{z) 



u 



2{N+k+l+l) 



(29) 



kJ=l 



where the choice of branch is determined such that F{z,{]) = 4:/{N — 1). Now equation 
( !28|) takes the form 



2k ,F{z;u)G{z;u,W) 



k=l 



(30) 
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where 

and the choice of the "+" sign in equation fl5Ul) determines the branch in the definition 
of u (essentially u^'^^ = y~^^^^^^'^). Now we introduce one more variable v such that the 
expansions for G mW and v agree to first order. That is, we define v such that 

r 4m2(A'+1)U/^ 2^2(^+1)^; 



giving 



2(Ar+l) 



Since ?/ ^ oo as 2; approaches along 7, it follows that u ^ along 7. From the 
boundedness of W and the definition of v ( |3Ti) . it follows that v is bounded on 7. We 
will now construct an initial value problem for z and v as functions of u that is regular 
at M = 0. 

Equation (!30|) takes the form 



2/ = -o >^ + o /JJ + 7^7 V^- 



2 ^ 2m^+i 

fc=i V ' y 

Recalling that y = u'"^, we have 

^ = n-(^-2)j(^;n,t;), (33) 

where J is analytic in its arguments at (2:0; 0, k), for any k G C, and J{zo] 0, k) = 1 — N. 

Now we will construct an equation for v'. Differentiating equation fl32l) with respect 
to 2; gives 

iy'= + v'. (34) 



Dividing equation ( l29i) by ^2(^+1) , differentiating with respect to z and using the reciprocal 
of equation ( l33l) to replace du/dz on the right side, shows that 



^t2(JV+l) 

F(z;m)2 



K{z;u,v), (35) 



where i^' is analytic in its arguments around {z; u, v) = {zq; 0, n), for any k E C Finally, 
from equation f|T8l) with S{z,y) = 0, we have 

/ „2(Af+l) \ 

W^' = -P(z, u')(v + — -^v' - 2Q{z, u^)u-^''+''>J{z- u, v) + R{z, u^) 

\ F[z]uY J 

^'^L{z;u,v), (36) 
16 



u 



where L is analytic in its arguments around (z;-u,f) = (2o;0,k), for any k G C. Using 
equations (135!) and (136|) in equation (!34l) shows that f ' is given by an equation of the form 



where H is analytic in its arguments around (z; m, f ) = (zq; 0, ft), for any « G C. Hence 

- = --=«J(.;«,.)i7(.;«,.). (37) 

Note that the system (|33ll and ( 1371) is regular at {z;u,v) = {zo]0,k). Since u —>■ 
and V is bounded on 7 near 2:0, it follows from Painleve's lemma (lemma [3]) that z and v 
are analytic functions of m in a neighbourhood of u = 0. Recall that J{zq; 0, = 1 — N. 
Hence from equation fl33l) . we see that z has a series expansion in u of the form 

00 

^ = ;,0 + ^«„w"+A^-l, 
n=0 

where ao = 1. Subtracting zq and taking the (A^ — l)st root gives 



00 



\ n=l 

for some choice of branch on the left side. Inverting this series shows that u is a power 
series in {z — l)i/(^~^) of the form 

00 



n=l 



where Aq = 1. It follows that y has a series expansion of the required form. 

Next we consider the case in which = 2K + 1 is odd. Returning to equation fl28l) . 
we let u = 1/y, giving 

k=l 

where e = ±1, F{z, u) is defined near {zq, 0) by 

. 2K , s 2K 

F{z; ur = ^^+8j2 ^^u'^-'^' + J2 b,{zMzy^^'^^^', (39) 



K"^ ^ k 

k=l k,l=l 



F{z, 0) = 2/K and 



17 



Note that the e = ±1 in equation (l38l) cannot be absorbed into a suitable choice of branch 
without affecting the value of F or G at m = 0. Following a similar argument as in the 
even N case, we define v by 

We then find that z and v satisfy a system of differential equations of the form 

- = u^-'Aiz;u,v), (40) 

^ = B{z;u,v), (41) 
au 

where A and -B are analytic in [z; u, v) at (2:0; 0, k), for any n in C, and A{z; 0, k) = eK. 
The existence of a solution of the appropriate form again follows from Painleve's lemma 
(lemma [3]). In particular, we see that z — zq ^ eu^ , so m ~ co{z — zqY^^ , where Cq = e. 

Case 3: A = 0. 

In this case the segments of the curve 7 on which y is small can be deformed through a 
region on which y is analytic such that y is bounded away from on the resulting curve. 
The proof is essentially identical to Shimomura ([26], lemma 2.2 and "case (iii)" following 
remark 3.1), see also [28] . 

Proof of theorem [2], part 2 The fact that these series converge follows from the exis- 
tence of the regular systems of equations for u and v (i.e., equations (|53H57|) in the even 
case and equations (jlOHlI]) in the odd N case). □ 



4 Singularity clustering along infinite length paths 

So far we have only considered singularities that can be reached by analytic continuation 
along finite length paths. The following proof is identical to that of the corresponding 
result in Smith [29] concerning equation (j4]) with deg/ > degg. 

Proof of theorem [5], part 4 

Let y be any solution of equation that can be analytically continued along a curve 
7 C up to but not including the point zq. Let D be an open disk in Q centred at zq. 
Without loss of generality, we assume that 7 C -D. Let zi be any point on 7. There exists 
a straight line segment / G D starting from zi such that analytic continuation of y along 
/ ends at a singularity Z2 of y in D. If this were not the case then y would be analytic 
throughout D, contradicting the assumption that y is singular at Zq. Since the singularity 
at Z2 can be reached by analytic continuation along a finite length curve /, it follows from 
part 3 of theorem [2] that it is of the algebraic type described. □ 

Next we summarize Smith's analysis [29] of equation (jl]). Although this example is 
not of the form of equation (jS]), it shows that the phenomenon described in part 4 of 
theorem [2] can indeed occur for some simple equations. 
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Smith showed that for any algebraic singularity of a solution of equation (jl]) at a point 
Zq G C, the function 

^z) := y\z) + viz)' 

is bounded in a neighbourhood of Zq. Furthermore, he showed that the only possible 
singularities of any solution of equation (jlj) are algebraic singularities or points at which 
such singularities accumulate along infinite length curves as described above. He shows 
that such an accumulation of algebraic singularities actually occurs by constructing a path 
7 ending at a finite point in the complex plane such that $ tends to infinity along 7. 

Smith found the following parametric representation of a one-parameter family of 
solutions of equation (jlj): 



viz) = 2-1/331/6^1/6 
. - ^-1/2 



-^-2/3(2; 



1/2 



r 









(42) 



. ^1/3(3;) . 
1/2 

(43) 



and 



X 



where Xq is an arbitrary constant and J,, is the Bessel function of order u. Smith showed 
that as X traces the path 

X = t — 1 cos 2 [t , zn < t < 00, 

V 6 4/ ' 

z traces a path 7 ending at a finite point z^: in the finite complex plane. Since ^{z) — >■ 00 
along this path, z^ must be an accumulation point of algebraic singularities. 

The kind of accumulation of singularities in Smith's example is far more complicated 
than the accumulation of poles in Painleve's example (j2]). The accumulation point of poles 
in Painleve's example can be reached by analytic continuation along finite length curves, 
so the standard methods used in the previous section show that no such accumulation 
can occur in many equations. 

Also, the more complicated accumulation of singularities along infinite length curves 
described by Smith's example (and allowed for in part 4 of theorem |2j) does not arise as a 
issue in the standard proofs that solutions of equations such as the first Painleve equation 
f fTTj) are meromorphic. In these proofs one assumes that a solution y is analytic at a point 
Zq and, arguing by contradiction, let a be the location of the nearest singularity that is 
not a pole. Since a can be reached by analytic continuation along a finite length curve, 
the standard arguments show that y is either regular or has a pole at a, i.e., we obtain a 
contradiction. It is only when the finite distance singularities are branch points that the 
accumulation in the finite plane along infinite length curves becomes a possibility. 

5 Discussion 

In this paper we have shown that when equation (jSj) admits "enough" formal algebraic 
series solutions at movable singularities, then such singularities are the only ones that can 



19 



be reached by analytic continuation along curves of finite length. The existence of these 
formal series solutions is equivalent to one {N even) or two (A^ odd) resonance conditions 
being satisfied. When these conditions are not satisfied, it is straightforward to show that 
there are formal series solutions of equation (|T3l) the form 



where Cq is a constant satisfying Cg ~ =1 and each c„ is a polynomial. If one tries to 
extend the proof above to include series of this form, one can construct a function W, 
which now includes a factor involving logy. Similarly, the corresponding u-v systems 
( !33] - [37l) and ( I401 - I4T1) have terms involving \ogu. In this setting, it is not sufficient to make 
estimates involving \y\ alone as one also needs to control the argument of y. 

Whether algebraic branch points in solutions of equation can accumulate in the 
manner described by part 4 for theorem [2] remains an open question. 

Acknowledgments 

The research reported in this paper was supported by a project grant from the Lever- 
hulme Trust. The second author was supported through and EPSRC Advanced Research 
Fellowship. We also acknowledge the support of the European Commission's Framework 
6 ENIGMA Network and the European Science Foundation's MISGAM Network. 

References 

[1] M.J. Ablowitz, A. Ramani, and H. Segur. Nonlinear evolution equations and ordinary 
differential equations of Painleve type. Lett. Nuovo Cim., 23:333-338, 1978. 

[2] M.J. Ablowitz, A. Ramani, and H. Segur. A connection between nonlinear evolution 
equations and ordinary differential equations of P-type. I and II. J. Math. Phys., 
21:715-721, 1006-1015, 1980. 

[3] M.J. Ablowitz and H. Segur. Exact linearization of a Painleve transcendent. Phys. 
Rev. Lett, 38:1103-1106, 1977. 

[4] J. Chazy. Sur les equations different ielles dont I'integrale generale est uniforme et 
admet des singularities essentielles mobiles. C.R. Acad. Sc. Paris, 149:563-565, 1909. 

[5] J. Chazy. Sur les equations differentielles dont I'integrale generale possede une 
coupure essentielle mobile. C.R. Acad. Sc. Paris, 150:456-458, 1910. 

[6] N. P. Erugin. Theory of moving singular points of second order equations. II. 
Differencial' nye Uravnenija, 12(4):579-598, 771, 1976. 

[7] A.S. Fokas, A.R. Its, A. A. Kapaev, and V.Yu. Novokshenov. Painleve transcendents, 
the Riemann-Hilbert approach, volume 128 of Mathematical Surveys and Monographs. 
American Mathematical Society, Providence, RI, 2006. 



oo 




n=0 



20 



L. Fuchs. Uber Different ialgleichungen deren Integrale Feste Verzweigungspunkte 
Besitzen. Sitz. Akad. Wzss. Berlin, 32:669-720, 1884. 

[9] V. V. Golubev. Lekcii po analiticeskoi teorii differencial! nyh uravnenii. Gosudarstv. 
Izdat. Tehn.-Teor. Lit., Moscow-Leningrad, 1950. 2d ed.]. 

[10] V.L Gromak, L Laine, and S. Shimomura. Painleve differential equations in the 
complex plane, de Gruyter, Berlin, 2002. 

[11] E. Hille. Ordinary differential equations in the complex domain. Wiley, New York, 
1976. 

[12] A. Hinkkanen and I. Laine. Solutions of the first and second Painleve equations are 
meromorphic. J. Anal. Math., 79:345-377, 1999. 

[13] M. Hukuhara. Meromorfeco de la solvo de (F) y" = Gy"^ + x. Lecture notes in 
Esperanto given at the University of Tokyo in 1960. Published by K. Okamoto and 
K. Takano in 2001 in ED. 



[14] E.L. Ince. Ordinary differential equations. Dover, New York, 1956. 

[15] N. Joshi and M.D. Kruskal. A direct proof that solutions of the six Painleve equations 
have no movable singularities except poles. Stud. Appl. Math., 93:187-207, 1994. 

[16] S. Kowalevski. Sur le probleme de la rotation d'un corps solid autour d'un point fixe. 
Acta Math., 12:177-232, 1889. 

[17] S. Kowalevski. Sur une propriete d'un systeme d'equations differentielles qui definit 
la rotation d'un corps solide autour d'un point fixe. Acta Math., 14:81-93, 1889. 

[18] B. Malgrange. Sur les deformations isomonodromique. I: Singularite regulieres. 
Seminaire de I'Ecole Morm. Sup. Birkhauser, Basel, Boston, Berlin, 1982. 

[19] B. M. McCoy, C. A. Tracy, and T. T. Wu. Painleve functions of the third kind. J. 
Mathematical Phys., 18(5):1058-1092, 1977. 

[20] T. Miwa. Painleve property of monodromy preserving deformation equations and 
the analyticity of r-functions. Puhl. Res. Inst. Math. Sci., 17:703-721, 1981. 

[21] K. Okamoto and K. Takano. The proof of the Painleve property by Masuo Hukuhara. 
Funkcial. Ekvac, 44(2):201-217, 2001. 

[22] P. Painleve. Sur les equations differentielles du premier ordre. C.R. Acad. Sc. Paris, 
107:221-224, 320-323, 724-726, 1888. 

[23] A. Ramani, B. Dorizzi, and B. Grammaticos. Painleve conjecture revisited. Phys. 
Rev. Lett, 49(21):1539-1541, 1982. 

[24] A. Ramani, B. Grammaticos, and T. Bountis. The Painleve property and singularity 
analysis of integrable and non-integrable systems. Phys. Reports, 180:159-245, 1989. 



21 



[25] M. Sato, T. Miwa, and M. Jimbo. Studies on holonomic quantum fields. VII. Proc. 
Japan Acad. Ser. A Math. Sci, 54(2):36-41, 1978. 

[26] S. Shimomura. Proofs of the Painleve property for all Painleve equations. Japan. J. 
Math, 29:159-180, 2003. 

[27] S. Shimomura. On second order nonlinear differential equations with the quasi- 
Painleve property II. RIMS Kokyuroku, 1424:177-183, 2005. 

[28] S. Shimomura. A class of differential equations of Pl-type with the quasi-Painleve 
property. Ann. Mat. Pura AppL, 186:267-280, 2007. 

[29] R.A. Smith. On the singularities in the complex plane of the solutions of 
y" + y'f{y) + giv) = P{x)- Proc. London Math. Soc, 3:498-512, 1953. 

[30] N. Steinmetz. On Painleve's equations I, II and IV. J. Anal. Math., 82:363-377, 
2000. 

[31] T. T. Wu, B. M. McCoy, C. A. Tracy, and E. Barouch. Spin-spin correlation functions 
for the two-dimensional Ising model: Exact theory in the scaling region. Phys. Rev. 
B, 13(l):316-374, 1976. 



22 



